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ABSTRACT 

We perform a calculation of the cross section for neutralino-neutralino annihi- 
lation into two photons and apply it to dark matter in the galactic halo to find the 
counting rate in a large gamma ray detector such as EGRET or ASTROGAM. 
Combining constraints from particle accelerators with the requirement that the 
neutralinos make up the dark matter we find that rates of over a few dozen events 
per year are unlikely. We list the assumptions that go into our conclusions and 
suggest other particle dark matter candidates which could give larger and perhaps 
observable signals. 
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I. Introduction 


It is possible that the dark matter (DM) which exists in the galactic halo 
consists of some, as yet undiscovered, elementary particle. If this is true and the 
particle is of the WIMP type (mass greater than 1 GeV and in thermal equi- 
librium at one time), then the particles may be detectable. 1 Schemes involving 
direct detection, through DM elastic collisions with nuclei, as well as indirect 
detection by observing DM particle-antiparticle annihilation products have been 
discussed and already limits have been placed on some popular DM particle 
candidates. 1 Since annihilation rates are proportional to the square of the DM 
number density, the strongest indirect limits have come from annihilation in the 
body of the Sun or Earth, where large density enhancements are likely, though in 
these cases it is only the neutrino annihilation products which can be observed. 1 . 

Detection of antiprotons, positrons and photons from annihilation in the 

galactic halo has been discussed 1 and found to be difficult; however, interest in 

2 

photons has been rekindled lately by the suggestion of Bergstrom and Snellman 
that annihilation into two photons could give a strong, extremely narrow line (at 
an energy equal to the mass of the DM particle) which would stand out against 
the diffuse gamma ray background. In fact, the relevant region of the gamma 
ray spectrum (> 1 GeV) will be measured for the first time in the very near 
future by the EGRET instrument on the GRO satellite 3 and larger devices with 
excellent angular and energy resolution such as ASTROGAM have recently been 
proposed. The mechanism of Bergstrom and Snellman is similar to the earlier 
suggestion of Srednicki, Theisen and Silk 5 that annihilation into charmonium 
plus photon would give an observable line, but is claimed to take place at a 
larger rate. 

Given the potential importance of the detection of such a line, we thought 
it worthwhile to check the Bergstrom and Snellman result, which was found 
using an approximation, by performing a more complete calculation. We have 
computed to leading order the gauge invariant set of one-loop box diagrams 
relevant for the process and find, for photinos or neutralinos, an annihilation 
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cross section which is more than an order of magnitude larger than the Bergstrom 
and Snellman result. In the pure photino limit this is in agreement with a very 
recent calculation of Rudaz, which was done following a different procedure. We 
also extend the result of Ref. 6 to the more realistic case of a general neutralino 
(generic combination of photino, zino and higgsino), taking into account the effect 
of the virtual exchange of possibly light Higgs bosons and charginos. Inclusion 
of these later particles can give further (but not large) enhancements. 

However, even with the enhanced cross section, the counting rate predicted for 
photino annihilation in a detector such as the GRO EGRET (or even the larger 
ASTROGAM) is found to quite small, probably less than a few dozen events 
per year and perhaps below the level of observability. The rate for the general 
neutralino is for most of the parameter space less than or equal to that for the pure 
photino special case. These conclusions come about, in part, because accelerator 
limits on photinos and squarks, in conjunction with limits from photino relic 
abundance, rule out the regions of parameter space where substantial annihilation 
would occur. Similar conclusions have been reached (though without the use of 
accelerator constraints) in a very recent paper by Bouquet, Salati, and Silk who 
used Rudaz’s cross section. 

It is important to keep in mind, however, that many different particle models 
can explain the dark matter and that some of these may evade accelerator bounds 
and perhaps produce observable signals. We give some examples of such models. 
In addition, uncertainties in astrophysical quantities such as the density and 
distribution of dark matter and uncertainties in the particle physics may make 
observable signals possible. 

The plan of this paper is as follows: In Section II we find the cross section 
by computing explicitly a gauge invariant subset of the contributing Feynman 
diagrams and expanding our result to leading order in 1/M* (Mj is the mass of 
the exchanged scalar fermion). The cross section for the generic neutralino (com* 
bination of photino, higgsino and zino) is given, though the large enhancement 
exists only for the photino/zino component. We then use the effective interaction 
technique to include the effects of chargino and Higgs exchange and discuss the 
limits of validity of our results. 
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In Section III we apply the cross section in the photino limit to find the flux 
in a gamma ray telescope such as ASTROGAM or EGRET. We pick values of 
the scalar fermion masses consistent with photino dark matter (.1 < Qj < 1), 
and show how this, along with with accelerator constraints from the Fermilab 
CDF experiment, the SLAC ASP experiment and TRISTAN experiments make 
rates of more than 40 events per year in the ASTROGAM detector (6 events/year 
for GRO) unlikely. We also display results for the general neutralino and show 
how these are in general dominated by the photino component. The Higgs and 
chargino contributions can be important in some regions of parameter space but 
in general no large enhancements are to be expected. The diffuse gamma ray 
background is also briefly discussed. 

In Section IV we consider models other than the minimal supersymmetric 
model. One simple model gives rise to a new dark matter candidate (the “luxino”) 
which would give easily observable fluxes of annihilation gamma rays even in the- 
EGRET detector, while avoiding accelerator constraints and being unobservable 
via other direct and indirect detection schemes. We point out that the “magnino” 
of Raby and West 8 would also produce an observable gamma ray line. 

In Section V we summarize the paper and make explicit various uncertainties 
and assumptions contained in our conclusions. We stress the importance of a 
high resolution gamma ray telescope in resolving some of these issues. Details of 
the cross section are displayed in the Appendix. 

II. Cross Section: XX — *■ 77 

In the minimal supersymmetric model which we consider there are four 
neutralinos which are linear combinations of the supersymmetric partners of 
the photon, Z ° , and two neutral Higgs bosons. The lightest of these makes 
an excellent candidate for the dark matter 10 ’ 11 ’ 12 and includes the pure photino 
and pure higgsino as special cases. We denote the lightest neutralino as x = 
ZuB + ZnW* + Z\zH^ + Z14H2 , where the axe the elements of the real or- 
thogonal matrix which diagonalizes the neutralino mass matrix. A pure photino 
has Z\\ = cos 0\VyZi2 = sin 9 w, Z13 = Z 14 = 0, while a pure higgsino can have 
Zu = Z 12 = 0, Z 13 = sin /?, and Z 14 = cos/3. 
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The neutralino masses and the are fully determined by four parame- 

ters: tan (3, /i, M, and M', where tan/3 = vi/v\ is the ratio of Higgs vacuum 
expectation values, M and M' are gaugino soft supersymmetry breaking mass 
parameters, and is a supersymmetric Higgs mass. Throughout, we make the 
standard simplification M' = | M tan 2 0w. The other relevant free parameters 
of the supersy mm etric model are the masses of the exchanged scalar fermions 
(squarks and sleptons). 

The annihilation into two photons takes place at the one-loop level through 
diagrams such those of in Fig. 1. (Additional diagrams with exchanged external 
photon and neutralino legs are not shown, for a total of fourteen diagrams.) 
Fig. 1 shows a gauge invariant subset consisting of scalar fermion plus fermion 
exchanges. Additional contributions from Z° plus fermion exchange are shown 
in Fig. 2, and additional contributions coming from chargino and Higgs exchange' 
will be discussed later. As will become apparent when we discuss rates in a 
gamma ray detector, we are interested in cases where the neutralino is less massive 
than the weak scale (m\y). We will therefore simplify the calculation by making 
an expansion in ( m x /mw ) 2 or ( m x /Mj ) 2 , where My is the sfermion mass. 

First consider the calculation of the amplitude given by the diagr ams of Fig. 1. 
We wrote the fourteen amplitudes for the general neutralino using the Feynman 
rides given in Refs. 9 and 13, and checked gauge invariance. We then performed 
the loop integrals after expressing them as three dimensional integrals over Feyn- 
man parameters in the standard way. At this point the expression for the total 
amplitude was quite lengthy and to simplify the calculation (and for the reason 
given above) we expanded in l/M? before evaluating the Feynman parameter 
integrals. The exception was the diagram of Fig. Id, where the integrand has 
a singular expansion in 1/My. Here we explicitly evaluated the parametric in- 
tegrals and then made the expansion, a rather lengthy calculation. Collecting 
terms and using several gamma matrix identities we found that all the gauge 
non-invariant pieces canceled and have a simple and gauge invariant expression 
for the amplitude 
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where pi, p 2 (fci, fo) are the four-momenta of the neutralinos (photons), ei and 
t 2 are the photon polarizations, m,j, Qi , and T 3i are the mass, charge and third 
component of weak isospin (±| for doublets and 0 for singlets) of the fermion in 
the loop, 3 is the Mandelstam variable, m x is the mass of the neutralino, and 


I(x) 
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( 2 ) 


Squaring this and taking the extreme non-relativistic limit valid here (the relative, 
velocity, v re j, is about 10“ 3 for galactic dark matter) we can write the resulting 
cross section for neutralino annihilation into two photons as 


O-Urel 


a 2 Gl 
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( 3 ) 


where Gf is the Fermi constant, a is the fine structure constant, ci is a color 
factor (3 for quarks, 1 for leptons), x{ = m x /rrii , and for the fermion plus sfermion 
exchanges considered in Fig. 1: 

4 =^f { Z n tan 0 w (Qi - T? L ) + Z n T? L ) 2 ^ 

b{ =0. 

It is important to note that the sum must be taken over both left and right 
chiral fermions. This allows the formula to be used with non- degenerate left and 
right sfermions. Also note that m x in eq. (3) carries the sign of the neutralino 
mass eigenvalue. 
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We have written the cross section, eq. (3), in this general form so that the 
contributions from Z° exchange as well as chargino and Higgs exchange can be 
included by just extending the sum to new objects. For example, the Z° exchange 
diagrams of Fig. 2 give a contribution 

A? =(Zl 4 

Bf = 0 , 


zh ) 2 
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m 2 z — 4 m x 2 
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— 4 m x 2 — 
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( 5 ) 


where Tz « 2.5 GeV is the width of the Z° and in this case the sum over the 
right chiral fermions gives zero. The contribution of the diagram of Fig. 2 where 
the fermion loop is replaced by a loop of W* bosons or squarks and sleptons will 
be neglected since it does not contribute to leading order in 1/m^- or 1/M?. 

In the pure photino limit, eq. (3) reduces to 


<rv ie i = — m. 


4- M). ^ X? - 

* /* ‘ 


( 6 ) 


(Note the sum is over both left and right sfermions, which gives a factor of 
four if they are degenerate.) A pure higgsino limit can also easily be obtained. 
Bergstrom, et al. 2 give a formula similar to eq. (6) with I(x)/x? instead of 
(1 — J(x)/x?). Eq. (6) is in agreement with the recent result of Rudaz, 6 and as 
shown by him is more than an order of magnitude larger than the earlier estimate 
of Bergstrom, et al.. 

Both the results of Ref. 2 and Ref. 6 were derived by using an effective 
interaction Lagrangian and computing the triangular fermion anomaly diagram. 
The disagreement comes from the way gauge invariance was imposed. Here, we 
took a slightly different and unambiguous approach. We performed the loop 
integrals of the box diagrams of Fig. 1, where no anomaly is present, explicitly 
retaining gauge invariance throughout. To leading order in the 1/M? expansion, 
the two procedures should give the same result, if performed correctly. In fact, 
only the diagram of Fig. Id (plus similar ones with exchanged external legs) can 
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give rise to a gauge invariant amplitude of order 1/M?. This can be understood 
by a dimensional argument and by noting that a gauge invariant amplitude must 
be at least quadratic in the external momenta, since the photons can appear in 
the effective interaction only through their field strength. Having shown now 
that the (much simpler) procedure of Rudaz leads to the correct answer, we can 
include the additional chargino diagrams using his method. 

In the minimal supersymmetric model there are two charginos, the super- 
symmetric partners of the and charged Higgs bosons, one of which is always 
heavier than the (and therefore contributes only to higher order in our ex- 
pansion). Their masses and mixing angles are determined by the parameters 
tan/?, M, and /*, defined earlier, and the lightest one can make an important 
contribution to our amplitude. Replacing the sfermion-fermion-neutralino ver- 
tex with the W r± -chargino-neutralino vertex (see Refs. 9 and 13 for the Feynman 
rules) we find a chargino plus W ± contribution given by the diagrams of Fig. 1, 
with the and chargino taking the places of the sfermions and fermions. As 
shown in the Appendix, the chargino plus W ± contribution adds another term 
to the sum in eq. (3): 

A **+) = “ 2 {0 2 l + Or) 


B^^-SOlOr, 


(7) 


where Ol = —ZuVu/y/2 + ^12^11> Or = Z\iU\i/^2 -f £ 12 ^ 11 , and the Uij and 
Vij are the matrices which diagonalize the chargino mass matrices (see Ref. 9 
appendices). When using eq. (7) in eq. (3), m,- = m x (+) and the color factor and 
charge are unity. 

A diagram similar to Fig. 2 with the fermions in the loop replaced by the 
lightest chargino likewise contributes: 


A #+) — g(^14 “ ^13)(^12 - ^12) f m 2 _ 

' z 


ra| — 4m x 2 


4 m x 2 — *r^m^ 
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( 8 ) 


>*<+) 


= 0 . 


Finally we consider the exchange of Higgs bosons. There are two scalar 
Higgs bosons, one pseudoscalar Higgs and two charged Higgs’ in the minimal 
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supersymmetric model. The scalar Higgs contribution vanishes in the extreme 
non-relativistic limit relevant here, and the charged Higgs’ contributions are small 
since they are heavier than the W ± . The pseudoscalar Higgs contributes via 
diagrams similar to those of Fig. 2 with the Z° replaced by the Higgs boson. For 
the pseudoscalar Higgs plus fermion exchange we have 

A? =0 

B? =2 2 mt T 2 (Zn tan 9 W - Z n )(Z u cos/3 - Z 13 sin/3)r<, (9) 

rrig — 

where rj = cot/3 for up type fermions and r{ = tan/3 for down type. We have 
included an extra factor of 2 in eq. (9) so that the sum x may be taken over the 
fermions only once (previously it had to be taken over both left and right chiral 
fermions). For the light chargino plus Higgs we also have a contribution: 

A? M =0 

2 

B$+) — 2 7^~j( Zl1 t&n &w - Zi2)(Zu cos/3 - Z i 3 sin/3)(l/n Vu cos /3 

771 ^ ^ 771 ^ 

+ U 12 V U sin/3). 

( 10 ) 

In the case of pseudoscalar Higgs boson exchange there are no loops of bosons 
or squarks and sleptons, so eqs. (10) and (9) give the complete result. 

It is important to note that our formulas are valid only in the limit m x 

Given, as shown in the next section, the low detection rates predicted 
for heavier neutralinos we feel that at this time the complete result would not be 
worth the effort and so we content ourselves with displaying results for m x < 40 
GeV. While heavier neutralinos may well make up the dark matter, it seems 
unlikely that they will be detectable via their gamma ray lines in the next few 
years. This concludes our discussion of the cross section. 

III. Rate in a Detector 

For a spherical isothermal galactic halo of core radius a, and a local dark 
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matter density the flux in a gamma ray detector can be written 




14,15 


d(j> 

dCl 


2 2 (< n, )ann a ( 1 + a2 ) 2 
Phal0 87rm x 2 03 


J(b,l ), 


( 11 ) 


where a = ro/a, ro is the distance of the Sun from the galactic center, (3 = 
(1 -f a 2 + a 2 cos 2 6 cos 2 Z) 1 / 2 , J(b,l) is an integral along the line of sight as a 
function of the galactic coordinates b and l, and we have included a factor of 
two for the two photons produced in each annihilation. We consider only high 
galactic latitudes (6 = 7r/2) where J « tt/2. 

Using 1 ro « 8.5 kpc, a ss 5.6 kpc, ph*lo 55=1 -3 GeV/cm 3 and measuring ( crv ) 
in units of 10 -26 cm 3 /sec, this becomes 


d<f> 110 (<rv) 26 events 
dfi (m^/GeV) 2 year sr cm 2 


( 12 ) 


For EGRET 3 with an effective area/field of view of 900 cm 2 sr this gives 


1.0 x 10 s (erv) 26 events 

l^GRO ~ "j , n \t\i > 

(m x /GeV) 2 year 

while ASTROGAM, 4 with an effective area/field of view of 7000 cm 2 sr gives a 
rate « 7.8 times larger. We shotild point out that there is great uncertainty in the 
astrophysical parameters a, r o and especially which enter into eq. (11), and 
that in addition, the halo is probably not an isothermal sphere. Taken together 
an order of magnitude uncertainty in the rate is probable. 

First consider pure photino dark matter. In order to evaluate eq. (13) using 
eq. (6) we need to choose values of the slepton and squark masses. These are 
free parameters of the minimal supersymmetric model most often considered, but 
fortunately the requirement that the relic abundance of photinos be consistent 
with the photino’s role as galactic dark matter constrains the masses considerably. 
If we require d.l < < 1, where f lx, is the ratio of the average photino density in 

the Universe to the critical density, and choose the value of the Hubble parameter 
h = 1/2, (50 km/sec/Mpc), only a fairly narrow range of sfermion masses is 
allowed. 
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In Fig. 3a we plot the values of the sfermion masses needed for Cljh 2 = 0.25 
(solid and short dashed lines), while in Fig. 3b we show the same for ft^h 2 = 
0.025. For the lines labeled “not split” we assumed that all squarks and sleptons 
(selectron, smuon, stau are the relevant sleptons) have a common mass M j. 
For the lines labeled “split” we took a common squark mass, Mq, and a common 
slepton mass, Me, but assumed Me = Mq/ 3. This is theoretically more attractive 
than the degenerate case. The slepton masses are labeled (e) (for selectron) and 
tire nearly the same as the common sfermion masses in the “not split” cases. This 
is to be expected since annihilation is dominated by the exchange of the lightest 
sfermion allowed. This will become important when we consider accelerator limits 
on the supersymmetric particles. 

The values of the sfermion masses, (and therefore the two-photon cross sec- 
tions) are found by requiring the indicated value of Clyh 2 for each value of the 
photino mass m,y. For each value of we solve the Boltzmann equation which' 
determines the relic abundance of photinos as they annihilated in the early uni- 
verse. We use an approximate, but very accurate (better than 5%) method which 
takes into account propagator momenta in the 77 -♦ qq cross sections as well as 
the changing degrees of freedom as the universe cools. See Refs. 17 and 11 for 
details of the method, cross sections and references. 

Using the sfermion masses from Fig. 3, we plot the counting rate from halo 
annihilation of pure photinos for the ASTROGAM detector in Fig. 4 and for 
GRO’s EGRET device in Fig. 5. Rates for both the split and degenerate cases 
are shown, although very little difference in the rates results. (Again, the rate 
is dominated by the lightest sfermions, i.e. the sleptons, and the Q® factor also 
favors the sleptons since they have unit charge.) 

For reference we also show a possible background rate as estimated by Stecker 
and Tylka . 18 The background rate is problematic since there are no measurements 
available at the relevant energies. Extrapolation from the 100 MeV region, where 
measurements have been made, seems unlikely to be accurate, and one therefore 
must rely upon theoretical calculations. The main source is thought to be photons 
from 7 T° decay , 18 where the pions are produced in cosmic ray interactions with 
the interstellar medium. A significant extragalactic component may (or may 
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not) exist and the background may be much lower in directions which have a low 
column density of ISM. We therefore consider the background estimates to be 
extremely uncertain. We simply integrate Stecker and Tylka’s estimate over the 
instrumental field of view and plot the number of events expected in an energy bin 
the size of the detector resolution (15% for EGRET , 3 and 1 % for ASTROGAM 4 
). The actual background could differ significantly from this when it is measured. 
See Refs. 7 and 18 for a more complete and careful discussion of background and 
how to extract a signal from it. 

More important at this point than the uncertain background is the low count- 
ing rates seen in Figs. 4 and 5. For the lowest value of Q^h 2 we consider (Fig. 4b), 
photinos of low mass do give substantial rates, hundreds of events per year being 
possible. Unfortunately, more careful examination of the squark and/or slep- 
ton masses required to produce this low Q^h 2 shows that photinos with masses 
less than around 11 GeV are inconsistent with either CDF or ASP experimental 
results. 

19 

This is shown in Fig. 3b where one long-dash line shows the CDF limit 
Mq < 74 GeV (90% c.l.) and the other long-dash line shows the ASP limit 20 
on the photino/selectron masses. The ASP line was found using cr(e + e~ — ► 
777 ) < 0.03 pb, and the relevant experimental acceptance. See Ref. 11 for the 
cross sections used, more details and further references. The TRISTAN limit 
(Mg > 26 GeV) does not further constrain the parameter space. 

Fig. 3b shows that for the case of degenerate squarks and sleptons (not split), 
the squark mass always falls below the CDF limit, ruling this case out. For the 
Mg = Mq / 3 (split) case, the squark masses evade the CDF limit for photinos 
above around 3 GeV in mass, but the selectrons fall below the ASP limits for 
photinos below around 11 GeV. Taken together we see that for the Ujh 2 = 
.025 case, photinos lighter than axound 11 GeV are ruled out, independent of 
squark/slepton splitting. Referring back to Fig. 4b we conclude that a rate larger 
than around-40 events per yeax in ASTROGAM is unlikely. We should also point 
out that for > 70 GeV, the squark/slepton masses necessary for Sl^h 2 = .025 
are less than m 5 .. This is inconsistent with our assumption of photino dark matter 
since only the lightest supersymmetric particle is expected to be stable. However, 
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as mentioned previously, tlie expansion used in calculation of the cross section 
breaks down before this and so we do not display these cases anyway. 

Fig. 3a shows that the limits from CDF and ASP for the = .25 case 
are much less stringent. For the degenerate case CDF allows photinos heavier 
than 9 or 10 GeV, while for the split case photinos heavier than 2 GeV are 
allowed by CDF. The combined CDF and ASP limits allow photinos heavier 
than around 5 GeV. However, Fig. 4a shows that the rates for the Cl^h 2 = .25 
case are correspondingly lower (for the same reason!), and again it seems that 
rates above a few events per year in ASTRO GAM are unlikely. 

As shown in Fig. 5, the rates expected in the soon to be launched GRO 
EGRET detector are substantially lower (and the background is higher), with 
less than 10 events per year expected if fl^h 2 = .025 and less than 1 per year if 
0 - r h 2 = .25. 

We should mention, at this point, that a way to evade the crucial ASP limit 
without lowering the photon event rate is to have a heavy selectron in conjunction 
with a light smuon or stau. The ASP cross section depends only on the selectron 
mass and therefore limits only this mass. However, while one may reasonably 
expect the colored squarks to split from the uncolored sleptons and even the stop 
squark to split from the other squarks due to the large top mass, there is no ready 
mechanism to split the sleptons among themselves, the electron, muon and tau 
have identical quantum numbers are are all relatively light, so we consider this 
possibility unlikely. 

Up to now we have considered only the pure photino special case of the 
neutralino, and only sfermion/fermion particles in the loops. This case is simple 
to display because the sfermion masses and rate are determined by the photino 
mass. For the generic neutralino there are more parameters (namely tan/3, fi 
and M, as discussed in Sec. II). However, we will now show that the simple case 
already displayed is very likely a rough upper limit to the rates possible for the 
generic neutralino. 

In Fig. 6 we show expected rates in the ASTROGAM detector for a set of 
supersymmetric parameters consistent with fijh 2 = .025. For comparison, we 
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plot (solid line) the rate for a pure photino from Fig. 4b (split case). All the 
points (x’s and boxes) are models with Mg = 120 GeV and Ms = 40 GeV, and 
therefore correspond in the pure photino limit to a single point on the solid line 
(near the big “blob”). The points are found by taking the above sfermion masses 
and values of tan (3 of 5, 2, 1.1, and .25, and then solving for all points (on a grid) 
in the /z, M plane which satisfy flj/i 2 = .025. Values of M and |/z| from 0 to 1 
TeV are considered, the x’s showing positive values of /z and the boxes negative 
values. 

To imagine the effect of including all of supersymmetric parameter space first 
fill in the area between all the points. This takes into account all values of tan/3, 
/z and M. Then “slide” the pure photino “blob” and attached filled area along 
the solid line, thus taking into account all possible values of squark masses. Note 
that the “blob” is the projection of the many points in supersymmetric parameter 
space which have large photino components. 

It is clear from Fig. 6 that the rate for the general neutralino seldom rises 
above the rate for a pure photino. This is to be expected since eq. (3) shows 
that any higgsino component is heavily suppressed with respect to any photino 
component, and the zino component is of the same order as the photino compo- 
nent. One might expect the chargino contribution (which is included in Fig. 6) 
to enhance the rate, since it can be very light. However, experimental searches at 
TRISTAN have ruled out charginos of mass less than 26 GeV and we therefore 
include only models (parameter values) which predict m x (+) > 26 GeV. Even 
apart from this, the chargino is really just another fermion in the loops of Figs. 1 
and 2, and since only one chargino is relevant, and it has no color factor, it cannot 
be expected to overwhelm the contributions from the many quarks and leptons. 
The small mass of the selectron (40 GeV) also enhances the effect of the sfermion 
exchange diagrams over the chargino exchange. In fact, for lower values of 
where larger sfermion masses are required, the role of the chargino is much more 
important (but the rates in a detector are smaller). It is interesting to note that 
one important effect of the chargino is to increase the rate for the pure higgsino 
case, which is usually thought to be quite suppressed. 

Next, we discuss the effects of the pseudoscalar Higgs boson. We have not 
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included the contributions from the Higgs in Fig. 6, even though it can be very 
light. A light Higgs might be thought to give a large enhancement due to the 
propagator {m\ — 4 m x 2 )~ 2 . However, the 4 m x 2 term dominates for light Higgs’, 
and except near the pole, the additional factor of ( m.i/mw ) 2 and the fact that 
the selectron mass is so low, makes the Higgs contribution small in comparison 
for the case under consideration. We did investigate numerically the effect of a 
light Higgs, finding that a substantial enhancement is present only close to the 
resonance. However, in these circumstances one expects neutralino annihilation 
in the early universe to be more efficient, diluting the relic density proportionally. 
Note that again, for lower the Higgs contribution can be dominant, but 

again the rates in a detector Me small. 

Finally, we point out that inclusion of the chargino and Higgs contributions, 
along with consideration of the general neutralino, make the very low rates ex- 
pected for the pure higgsino not nearly as likely. Detection of the general neu-. 
tralino is most likely easier than detection of a pure higgsino. Due to the low 
rates, however, a full exploration of the supersymmetric parameter space will not 
be conducted at this time; rather we will explore other models which have the 
potential for giving large and probably observable signals. 

IV. Models with Large Rates 

In this section we attempt to find particle models which evade existing particle 
accelerator constraints and give large gamma ray signals. The counting rate in a 
gamma ray detector from DM annihilation is low primarily for two reasons. First, 
the number density of DM particles in the halo is low, and the annihilation rate is 
proportional to the density squared. Second, the cross section for annihilation of 
neutral particles into two photons typically occurs only at the loop level, implying 
that the annihilation cross section is low. Now the total annihilation cross section 
of a candidate DM particle can be determined (at “freeze-out”) by requiring a 
cosmologically relevant relic abundance 

(<™>«n«~ nh 2 sec > ( 14 ) 

where some dependence upon the “p” vs. “s” wave annihilation has been left 
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out. The problem with the neutralino (and massive neutrino, etc.) is that the 
annihilation cross section which eq. (14) determines is dominated by the quark- 
antiquark and lepton-antilepton channels, leaving the two photon cross section 
smaller by factors of a 2 / 7r 2 . As a way around this, we consider a class of particles 
(called “luzinos ” , for “lux”: light) whose main annihilation channel in the early 
universe is into photons. Assuming this, eqs. (13) and (14) predict 


-Rgro ~ 


10 5 (mj/GeV) 2 events 
Cllh 2 year ’ 


(15) 


and a rate about eight times larger for ASTROGAM. Here mj and f 1/ are the 
luxino mass and relic energy. This substantial rate would be easily visible above 
background by the EGRET detector. 

As an existence proof that such particles are possible, consider a model with 
a neutral Majorana fermion luxino Z°, a pair of charged fermions and a neutral 
scalar <f>. Assuming the 1° is a singlet under the standard model gauge group, we 
introduce the interaction Lagrangian 


£, = Z°(a + &7s)/V + rV + + h.c.. (16) 


For simplicity, we will consider the case where and (f) have a common mass M 
(larger than the luxino mass mi). The luxino annihilates through the diagrams 
of Fig. 7 leading in the non-relativistic limit to the cross section 




77)v ~ 


4 <» 2m f t2 
7 r 3 M 6 


(-a' 2 


+ b ,2 ), 


(17) 


where we have taken M » mj. The other channels, l°l° — > (7 *,Z*) —* ff are 
not able to compete with l°l° —* 77, since 1° develops non- vanishing electroweak 
form factors only at the two loop level. This also implies that the luxino has a 
very feeble electroweak coupling with ordinary matter, which makes it invisible 
in direct dark matter detection experiments via nuclear recoil. If we take the 
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coupling constants 6 , a' and b' to be of order unity, the present luxino density is 
given roughly by 


Qlh 2 as 8 x 1CT 5 


(Af/GeV ) 6 

(mj/GeV ) 4 


(18) 


Dark matter luxinos in the mass range of 10 to 50 GeV require M in the range 
25 to 70 GeV. In order to make the model complete, we need to add to the 
Lagrangian, eq. (16) some interactions able to mediate the decay of the charged 
/ ± particles. One possibility is to introduce heavy (> TeV) charged scalars or 
gauge bosons which couple a ft-l 0 pair to ordinary fermions (/) and allow the 
process — ► l°ff. 

This example illustrates a particle which predicts a photon counting rate as 
large as eq. (15), but which is invisible in other accelerator and dark matter 
search schemes. 

Finally, we want to mention that large monochromatic gamma ray production 
can be expected in the model of the magnino, proposed in Ref. 8 to solve the solar 
neutrino and dark matter problems. In this model, the spin \ neutral magnino 
(m°) has a Yukawa interaction with a charged fermion m + and a Higgs scalar 
h + . The masses of m + and h + are almost degenerate with the magnino mass in 
order to provide it with a large magnetic moment and still escape experimental 
detection. The annihilation m°m —* 77 occurs through the diagrams of Fig. 1 
(substituting m + and h + in the internal lines). In the non- relativistic limit, we 
obtain a cross section of order 


a 2 a 4 

a(m°m° -» 77 ) 1 ; ~ g - 3 ^ 2 , (19) 

where g is the Yukawa coupling of order unity, and mo is the magnino mass, which 
must be in the range 5 to 10 GeV to solve the solar neutrino problem. Taking into 
account that, in order to solve the solar neutrino problem, the number density 
of anti-magninos must be less than the number density of magninos by a factor 
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of two or so, we obtain a maximum counting rate for the EGRET detector of 


-Rgro 


s, , „ , r ._, events 

5 x 10 5 (mo/5 GeV) . 

year 


(20) 
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It is important to recall, however, that a recent reanalysis of the e + e~ annihila- 
tion data from the Mark II detector at SLAC has left only a very narrow allowed 
region for the magnino mass parameters and the model seems disfavored. 

V. Conclusions and Discussion 


We have performed a careful calculation of the cross section for the annihila- 
tion of galactic neutralinos into two photons and shown that the rate of annihila- 
tion can be substantial, mainly due to the gaugino (photino and zino) component, 
dominating the rate into charmonium plus photon. Two ’ calculations of the 
this cross section differ by more than an order of magnitude and we have shown 


using a different procedure that the enhanced results of Rudaz are correct to first 
order in Mj 2 . We have included the effects of charginos and Higgs bosons in the 
loop diagrams and shown that while they can make important contributions, no 
further large enhancements are to be expected. 


Even with this greatly enhanced cross section, however, we find for canonical 
values of the galactic halo parameters, that less than a few dozen events per year 
will be seen in a gamma ray detector the size of ASTROGAM. This is due, in 
laxge part, to the constraints put on dark matter photinos by accelerator results 
from the CDF and ASP experiments. 


To see if other dark matter candidates might produce more interesting gamma 

ray signals we then introduced a new class of models (the “luxinos”) whose main 

annihil ation ch ann el is two photons. These can give much larger, probably ob- 

8 

servable signals. We pointed out that the magnino of Raby and West can also 
provide large signals. 

Since the data available on the identity of the dark matter is so scarce, it 
is important that high resolution gamma ray searches take place. These will 
complement the direct detection experiments and the search for annihilation 
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neutrinos in proton decay detectors. While we did not emphasize background, 

we should note that the intrinsic width of the gamma ray line is ~ 10 -3 and 

that a detector with an energy resolution of 10 -3 could cut the background by 

a factor of 100 over ASTROGAM. It is also the case that uncertainties in the 

distribution and density of dark matter could increase the rate substantially over 

our estimates. This is especially true if there is a large concentration of dark 

23 

matter near the galactic center. 

Finally we should note that our neutralino results are only valid to first order 
in m x IMj and that above around 40 GeV they we not to be trusted. While it 
is certainly possible that neutralinos more massive than this make up the dark 
matter, the photon counting rate is dropping rapidly and the substantial increase 
in cross section that would be needed to allow easy observation of these heavy 
neutralinos seems unlikely. 

In conclusion, we feel that while the outlook for neutralino detection in the 
near future via halo annihilation is rather pessimistic, the uncertainties in the 
predictions are large enough to warrant a careful search for sharp lines in the 
diffuse gamma ray background. 
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Appendix 

In this appendix we derive the Ai and Bi functions of eq. (3) using the effective 
interaction technique of ref. 6. The effective Lagrangian, valid for M j,myy >> m x , 
for the interaction between neutralinos and fermions (quarks or leptons) is 


C *ff = v/ 2 G FX7 / *75X + <*RPR)f, 


(Al) 
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where 


oil —(^i 4 — ^13 ){Qf sin 2 


T} l ) - 


m 


W 


M h 


Z\\ tan Ow{Qf ~ T} L ) + ZuT} 


i3£ 


1 2 


2 m? 

(xr -{Z\ A - Z\i)Qf sin 2 Ow + (Zu tan 8wQf) 2 , 

M fR 

(A2) 

Pi = (1 — 75)/2, Pr = (1 4- 75)/2, and all other symbols were defined in Sec. II. 
For the neutralino-chargino interaction we have 



= V2G f 


X7 m 75X X (+) 7 h(PlPl + 0rPr)x^ +) 

+ 05 (xx X (+) X (+) + X75X X (+) 75X (+) ) j ) 


(A3) 


where 

0L ={Zi 4 — Z i 3 ) ^—1 + sin 2 9\y + -Vj 2 2 ^ + 20£ 

Pr =(Z 2 u - ^ 13 ) ("1 + sin 2 0 W + - 20r 

05 =40iOr , 


(A4) 


and Ol and Or were defined just after eq. (7). 

The coupling of the neutralino current to the quark/lepton and chargino 
currents leads to an effective interaction between the two neutralinos and a pair 
of photons through the diagrams of Figs. 8. For a fermionic current /A/, where 
A = 1, 75, 7 a, 7a 75 for scalar, pseudoscalar, vector and axial vector couplings 
respectively, the effective Feynman rules for the triangular loop diagrams of Fig. 8 
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are 


scalar =* ^ Q}cimi[l + (1 - ^)/(*j)] 

pseudoscalar =e tiupa ' k\pk 2v t\ p C2<T — Q?CimjJ(x;) 

7T 

t 

vector =0 

axial vector = - e* u ' p<T ki^i/Clp^C^i + &2 )a~ Y'' <?iC;[l - 

7T Xj 


( 45 ) 


where /ki f 2 and ei f 2 are the momenta and polarization of the two photons, with 
{k\k 2 ) the dot product given by a metric with signature (1, —1, —1, —1). The sum 
runs over the fermions exchanged in the loop and here x? = (fci + k 2 ) 2 /(4m?). All 
other symbols were defined in Sec. II. Prom the above effective interactions and 
Feynman rules one can easily derive the cross sections to leading order in l/Mj 
or 1 /rriyy as presented in Sec. II. In the case of the loop of and charginos, the 
validity of the effective interaction procedure can be checked by working in the 
unitary gauge and retaining only the leading terms of the the W ± propagator. 
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Figure Captions 


1 . Feynman diagrams contributing to neutralino (x) annihilation into two 
photons ( 7 ). Additional diagrams with exchanged neutralino and photon 
external legs are not shown. The symbol / stands for either a left or right 
chiral sfermion (squark or slepton) and / for a left or right chiral fermion 
(quark or lepton). 

2. Feynman diagram for Z° exchange contribution to neutralino annihilation 
into two photons. The diagram with exchanged photon legs is not shown. 

3. Constraints on sfermion (squark and slepton) masses from relic abundance 
and particle accelerator experiments. As a function of the photino mass 
(m 7 ), the masses of the squarks (solid lines) and selectrons (short-dash 
line) required to give the indicated value of ft^h 2 are shown. Fig. 3a is for 
fl 7 h 2 = .25 and Fig. 3b for Q^h 2 = .025. Both the split (Mg = M$/Z) and 
not split (Mg = M§) cases are shown. One long-dash line shows the CDF 
lower limit on squark masses, while the other shows the ASP limit on the 
selectron mass. Regions below these lines are ruled out. 

4. Predicted counting rate of photons in the ASTROGAM detector from anni- 
hilation of galactic photinos as a function of the photino mass. The sfermion 
masses were found by requiring f hjh 2 = .25 (4a) or Q^h 2 — .025 (4b) (see 
Fig. 3). Both the split (Mg = M^j 3, short-dash line) and the not split 
(Mg = M$, solid line) cases are shown, as is a background counting rate 
estimate of Stecker and Tylka (long-dash). See Sec. Ill for more details. 



5. Same as fig. 4 for the GRO EGRET device. 


6. Photon counting rate for ASTROGAM detector from general neutralino 
annihilation. Chargino exchange is included and till points satisfy tl^h 2 — 
.025. See Sec. Ill for more details. 

7. Feynman diagram contributing to luxino (1°) annihilation into two photons. 
The diagram with exchanged photons is not shown. 

8. Feynman diagram for effective interaction of neutralinos and photons. * 
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